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Abstract 

Plumbagin (5-hydroxy-2-methyl-1,4-naphthoquinone) is a naphthoquinone derivative from the roots of plant Plumbago 
zeylanica and belongs to one of the largest and diverse groups of plant metabolites. The anticancer and antiproliferative 
activities of plumbagin have been observed in animal models as well as in cell cultures. Plumbagin exerts inhibitory effects 
on multiple cancer-signaling proteins, however, the binding mode and the molecular interactions have not yet been 
elucidated for most of these protein targets. The present study is the first attempt to provide structural insights into the 
binding mode of plumbagin to five cancer signaling proteins viz. PI3Ky, AKTl/PKBot, Bcl-2, NF-kB, and Stat3 using molecular 
docking and (un)binding simulation analysis. We validated plumbagin docking to these targets with previously known 
important residues. The study also identified and characterized various novel interacting residues of these targets which 
mediate the binding of plumbagin. Moreover, the exact modes of inhibition when multiple mode of inhibition existed was 
also shown. Results indicated that the engaging of these important interacting residues in plumbagin binding leads to 
inhibition of these cancer-signaling proteins which are key players in the pathogenesis of cancer and thereby ceases the 
progression of the disease. 
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Introduction 

Cancer is the leading cause of disease worldwide accounting for 
12.7 million new cases every year and this number is expected to 
rise to a whopping 26 million by 2030 [1]. Considering the impact 
on human health and economics, cancer presents a major 
challenge to the scientific world and there is a necessity to 
discover novel agents for the treatment of this disease. Several 
studies have been focused on naturally occurring chemical 
compounds which are known to possess cytotoxic effects and 
have the potential for killing cancer cells [2-3] . Plumbagin (PL; 5- 
hydroxy-2-methyl-l,4-naphthoquinone, Figure 1) is one such 
important compound, which is a naphthoquinone derivative 
identified from the roots of plant Plumbago zeylanica and belongs 
to one of the largest and diverse groups of plant metabolites [4—6]. 
The anticancer and antiproliferative activities of PL have been 
observed in animal models as well as in cell cultures [7-12]. The 
anticancer activities of PL have been shown against a wide variety 
of cancers including breast cancer [13], lung cancer [14—15], 
ovarian cancer [16], acute promyelocytic leukemia [17], and 
prostate cancer [18-19]. The inhibitory effect of PL was shown by 
alterations of various signaling pathways which play a crucial role 
in cancer cell proliferation, survival, invasion, and metastasis 



[14,16,18,20-22] through suppression of major signaling mole- 
cules such as nuclear factor-kappaB (NF-kB) [23], AKT/mTOR 
[13], and signal transducer and activator of transcription 3 (Stat3) 
[24]. 

In the current study, we considered five key signaling molecules 
namely phosphatidyl inositol-4,5-bisphosphate 3-kinase (PI3K), 
AKT (also known as PKB, protein kinase B), anti-apoptotic 
protein B cell lymphoma-2 (Bcl-2), NF-kB, and Stat3 which play 
important role in cancer development and progression. These 
molecules have already been reported [13,18,23,25] to show direct 
inhibition by PL. Owing to their key role in cancer development, 
these regulatory molecules are frequently targeted in anticancer 
therapy. 

The PI3K and AKT pathway, commonly known as PI3K/ 
AKT/mTOR pathway, is perhaps the most frequendy dysregu- 
lated pathway in human cancers [26—29]. Consequently, signifi- 
cant efforts have been made to generate inhibitors of the key 
kinases in PI3K pathway including PI3K, AKT, and mTOR in 
recent years [30] . Plumbagin targets PI3K in human breast cancer 
cells and drastically decreases the level of PI3K subunit p85 
causing downstream AKT/ mTOR pathway inhibition leading to 
growth arrest and cell death [13,25]. The PL target, AKT, plays 
an important role as antiapoptotic protein [31,32] and is a widely 
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Plumbagin 
(5-hydroxy-2-methyl-1,4-naphthoquinone) 

Figure 1. Two dimensional structure of plumbagin. The lUPAC 
name of the plumbagin is '5-hydroxy-2-methyl-l,4-naphthoquinone'. 
The carbon atoms are shown with the lUPAC numbering. The numbers 
are referred in the text wherever a particular atom of plumbagin is 
involved in molecular interaction with target proteins. 
doi:1 0.1 371/journal.pone.0087309.g001 

exploited for anticancer therapy [33-35]. Plumbagin inhibited 
AKT kinase activity and induced autophagy in two cancer cell 
lines, MDA-MB-231 and MCF-7 [13]. 

The signaling molecule NF-kB transcription factor plays a 
major role in the development and progression of various types of 
cancer [36-38]. A constitutive and continuous NF-kB activity is 
observed in various cancer tumors including lymphoid or myeloid 
tumors [39-42]. In addition, tumor cells often exploit NF-kB to 
achieve resistance to chemo- and radio-therapy [43]. Anti-cancer 
therapies have been developed by exploiting the inhibitors of NF- 
kB which may reduce the progression and development of the 
disease and/ or may improve the efficacy of conventional therapies 
[44]. In this regard, certain inhibitors were discovered which 
suppress activity of NF-kB by directly blocking its binding to DNA 
[45-47] . Using the similar mechanism, PL inhibits NF-kB activity 
by interfering with the binding of DNA to NF-kB as demonstrated 
by gel shift assay of nuclear proteins [23,25]. 

The Bcl-2 family of proteins are the principal regulators of 
programmed cell death or apoptosis [48,49]. The over-expression 
of Bcl-2 has been observed in various cancers which contributes to 
drug resistance [50] and enhances in-vivo cell survival [5 1] while 
reduction in Bcl-2 expression increases drug sensitivity [52]. 
Several reports [53-58] have explored the development of 



inhibitors of Bcl-2 protein as potential anti-cancer drugs. In this 
regard, plumbagin has also been shown to induce apoptosis by 
down-regulation and inactivation of Bcl-2 in human breast cancer 
ceUs [25]. 

The signaling molecule StatS is constitutively activated in many 
human cancers and has been widely exploited as therapeutic target 
for cancer therapy [59-63]. In this regards, plumbagin has been 
shown [18] to inhibit the DNA-binding activity of Stat3 in prostate 
cancer ceU lines, DU145, PC-3, and CWR22rvl. 

Although PL mediated inhibition of the aforementioned five 
signaling molecules has been reported as reviewed above, the 
binding mechanisms and the molecular interactions of PL with 
these key cancer regulatory molecules is apparently not known. 
Plumbagin docking has been reported previously with few protein 
molecules [16,64-68], however, the present study is the first 
attempt to investigate the structural and molecular details of PL 
binding against five key cancer signaling targets viz. PISKy, 
AKTl/PKBa, Bcl-2, NF-kB, and Stat3 using molecular docking 
and ligand (un)binding simulation approach. In addition, we also 
identified and ranked some important interacting residues of these 
target proteins. 

Materials and Methods 

Data retrieval 

The molecular structure of PL was retrieved from PubChem 
compound database with CID 10205. The 3-D structures of five 
chosen cancer signaling targets were retrieved from Protein Data 
Bank (PDB) and are: AKTl (PDB ID: 3096), Bcl-2(PDB ID: 
2021), NF-kB (PDB ID: 3GUT), PI3Ky (PDB ID: 3L54), and 
Stat3 (PDB ID: IBGl). The bound hgand was used as probe for 
the binding site grid generation. In case of StatS protein, the 
available structure is from mouse and differs from the human 
homolog in an amino acid residue at 760 position (Asp instead of 
Glu in human). The Stat3 structure contains 3D atomic 
coordinates until the amino acid residue at 716 position. 
Therefore, being 100% identical within the whole range of amino 
acids for which 3D coordinates are available, it will presumably 
serve as Stat3 structure for human as well. All the software tools 
used in this study are summarized in Table 1. 

Molecular docking 

Dock v.6.5 from the University of California, San Francisco was 
used for all computer simulations and docking of PL into the active 
site pockets of five chosen cancer drug targets was achieved [69] . 
The best conformation search strategy exploited in the present 
work was random conformation search, which utilized the existing 
Coulombic and Lennard-Jones grid-based scoring function. The 
whole process of docking involves multiple steps. Briefly, once the 



Table 1. The computational tools which were used in this study are presented with their source and application. 





Computational tool 


Availability 


Application 


Dock v.6.5 


http://dock.co m p b io . u csf.edu/ 


Molecular docking of ligand into the binding site of the 
receptor. 


PyMOL v.1.3 


http://www.pymol.org/ 


Molecular visualization and analysis 


LigPlot+ 


http://www.ebi.ac.uk/thornton-srv/software/LigPlus/ 


Molecular interaction between receptor and ligand. 


Naccess v.2.1.1 


http://www.bioinf.manchester.ac.uk/naccess/ 


Loss in accessible surface area of residues 


X-Score v.1.2.11 


http://sw16.im.med.umich.edu/software/xtool/ 


Binding energy and dissociation constant of ligand 


MoMA-LigPath 


http://moma.laas.fr 


Simulate protein-ligand (un)binding 
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Figure 2. Plumbagin (PL) binding to the cavity of PISKy. Panels A-F: show the (un)blndmg simulation phases of PL, 'A' is farthest from the 
binding site, 'E' is the closest to the binding site, and 'F' is the binding site phase. The hydrogen bonds are shown as green-dashed lines with 
indicated bond length and the residues involved in hydrophobic interactions are shown as red arcs. Those residues which are common to the last 
phase (F) are encircled. Panel G: Another representation for phase F. The whole protein is displayed in cartoon representation and the ligand 
molecules are in sticks; PL colored as green and the bound known inhibitor in blue. The interacting residues are labeled and shown as surface in 
different colors. 

doi:1 0.1 371 /journal.pone.0087309.g002 

potential site of interest is identified on the receptor, a grid is interest is calculated using a set of matching atom-sphere pairs, 
generated within the site and each grid point is considered as a To evaluate the orientation of the ligand within the site, a 
sphere center. The orientation of the ligand within the site of shape scoring function and/or a function approximating the 
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ligand-protein binding energy is used. The ligand-protein binding 
energy is an approximate sum of tlie van der Waals attractive, van 
der Waals dispersive, and Coulombic electrostatic energies. As a 
final step, to minimize the energy score, the orientation of the 
ligand is varied slighdy and evaluated. After the initial orientation 
and evaluation (scoring) of the ligand, a simplex minimization is 
used to reach the nearest local minimum of the energy score. The 
orientation of the ligand corresponding to the nearest local 
minimum energy score serves as the final docked ligand. 

Analysis of docked protein-ligand complex 

PyMOL V. 1.3 [70] was used to analyze and generate an 
illustration of whole protein-ligand complex. LigPlotH- v. 1.4. 3 
program [71-72] was used for analyzing the interaction of docked 
protein-ligand complex to check the polar and hydrophobic 
interactions between the receptor and ligand and illustrations of 
molecular interactions between PL and the chosen proteins were 
generated. To confirm the involvement of interacting residues 
obtained from LigPlot+, loss in Accessible Surface Area (ASA) was 
calculated. If a residue lost more than 10 ASA in the direction 
from the unbound to the bound state, it was considered to be 
involved in interaction [73]. The ASA of unbound protein and tiie 
protein-ligand complex were calculated using Naccess v.2. 1 . 1 [74] . 
The change in ASA (AASA) of the residue in the direction from 
unbound to bound state was calculated using the expression: 



AASA, = ASAf°"='" - ASA™*""-"^''"'' 

In addition to the Dock score (Grid score) obtained from Dock 
[69], the binding energies and dissociation constants were also 
calculated by X-Score v. 1.2.1 1 [75-76]. 

Protein-ligand (un)binding simulation 

MoIvIA-LigPath (http://moma.laas.fr), a web server which 
works on Molecular Motion Algorithms (MoMA) [77-78], 
simulates the ligand unbinding from the binding site of the 
protein to the surface of the protein. It considers the flexibility for 
the protein side-chains and the ligand, and involves only geometric 
constraints. This method provides mechanistic information about 
how ligand is driven to the binding site from the surface of the 
protein or from the binding site to the surface. It also provides 
snapshots of molecular interactions leading the ligand from surface 
of the protein to the binding site. In the process, it identifies those 
residues which despite being away from the binding site still play 
important role in ligand binding or in driving the ligand to the 
binding site. The docked molecular complex was subjected to 
(un)binding simulation using MoMA-LigPath. 

Results and Discussion 

Plumbagin docking and (un)binding simulation study of 
PI3Ky 

The snapshots of molecular interactions leading the ligand from 
surface of the protein to the active site are shown in Figure 2A-G. 
Final binding in the active site and the molecular interactions 
holding PL in the active site are shown (Figure 2F-G, Table 2). 
Plumbagin packs against the residues Ile-963, Met-953, Val-882, 
Ile-881, Glu-880, Ile-879, Tyr-867, Ile-831, Trp-812 and Met-804 
(Table 3). In the final docking phase of PI3K")', a hydrogen bond 
was formed between O-atom of 1-carbonyl group of PL and ot- 
amino group of Val-882 and other 34 hydrophobic interactions 
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Table 3. The plumbagin interacting residues of the five cancer signaling proteins are listed with the number of non 
contacts and the loss in Accessible Surface Area (ASA). 


-bonding 




Target 


Plumbagin Interacting residues 


No. of Non-bonded contacts 


AASA (in A^) 




PI3Ky (pi 10) 


Met-804 


1 


15.69^ 






Trp-812* 


2 


20.37* 






lle-831 


3 


24.31^ 






Tyr-867 


8 


6.19'' 






lle-879 


1 


12.18' 






Glu-880 


2 


4.92'° 






lle-881 


2 


7.61" 






Val-882 (H-bonding)* 


4 


12.45"^ 






Met-953» 


7 


31.39' 






lle-963 


4 


22.66^ 




AKTl 


Trp-80* 


20 


59.53' 






Ser-205 


1 


16.08^ 






Leu-210 


1 


19.41=^ 






Thr-21 1 


2 


12.79*^ 






Tyr-272 


1 


18.42^ 






Asp-292 


2 


18.36'' 




Bcl-2 


Thr-93 


1 


17.78^ 






Gln-96 


1 


18.31^ 






Ala-97 


2 


14.31^ 






Tyr-199* 


17 


47.08' 






Pro-201 


1 


16.46* 




NF-kB (p65) 


Tyr-36* 


4 


26.13^ 






Val-121 


1 


4.21^ 






Lys-122 


2 


7.40" 






Lys-123 (H-bonding)* 


4 


36.19' 






Asn-155 (H-bonding) 


0 


17.53^ 




Stat3 


Arg-382 


3 


17.63^ 






Glu-415 


1 


1 7.80^ 






Arg-417 (H-bonding)* 


2 


42.93' 






Cys-418* 


3 


1 7.74^ 






Gly^l 9* 


1 


6.91** 




The residues involved in hydrogen-bond formation are indicated by 'H-bonding' 

accessibility indicated by superscripts with the value of AASA. 

*The residues which were common across all phases of {un)binding simulation. 


in parentheses. The residues were also ranked on the basis of loss in 


solvent 



doi:l 0.1 371 /journal.pone.0087309.t003 



were from 10 different residues (Figure 2F, Table 2, Table 3). The 
Dock score was highly negative and binding energy was also 
comparable to that of the known inhibitor while the dissociation 
constant was approximately one tenth of the known inhibitor 
(Table 2). The residues Trp-812, Val-882, and Met-953 were 
common in all the phases of PL (un)binding simulation, showing 
their importance in initial recruitment of PL from surface of the 
protein to the deep binding site (Figure 2). Previously, the residues 
of PISKa subunit viz. Ile-800, Ile-848, Val-850, and Val-851, (the 
equivalent residues in PISKy are Ile-831, Ile-879, Ile-881 and Val- 
882 respectively) were also reported to play a role in binding of the 
inhibitor wortmannin. Specifically, the residue Val-85 1 (equivalent 
to Val-882 of PISKy) was involved in hydrogen bonding to 
wortmannin through its amide group [79]. Also in an earlier 
report [80], a highly potent inhibitor GSK2 126458 of PISKy was 
co-crystallized with the protein (PDB IDs: 3L54, 3L08) in which 



Val-882 was involved in hydrogen bonding. In another study [81], 
Val-882 through its backbone was also found to be involved in 
hydrogen bonding with H-bond acceptor in five different 
inhibitors i.e. wortmannin, LY294002, quercetin, myricetin, and 
staurosporine. The wortmannin was found to be packing itself 
against N-terminal lobe residues of PISKy viz. Ile-831, Ile-879, Ile- 
881, and Val-882 [81] which are overlapping with the PL binding 
site. The hydrogen-bonding feature by the residue Val-882 was 
reported to be conserved in all kinase-inhibitor complexes [82]. In 
this regard, our docking simulations of PL and PI3Ky were in 
agreement with the above studies e.g., the conserved H-bonding 
feature of Val-882 and the other common residues involved in PL 
binding. These consistent findings provided validation for our 
molecular docking and also showed that PL inhibits PI3Ky with 
the mechanism similar to those of other reported inhibitors. 
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Figure 3. Plumbagin (PL) binding to tlie cavity of AKT1. Panels A-F: show the (un)binding simulation phases of PL, 'A' is farthest from the 
binding site, 'E' is the closest to the binding site, and 'F' is the binding site phase. The residues involved in hydrophobic interactions are shown as red 
arcs. Those residues which are common to the last phase (F) are encircled. Panel G: Another representation for phase F. The whole protein is 
displayed in cartoon representation and the ligand molecules are in sticks; PL colored as green and the bound known inhibitor in blue. The 
interacting residues are labeled and shown as surface in different colors. 
doi:1 0.1 371/journal.pone.0087309.g003 
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Plumbagin docking and (un)binding simulation study of 
AKT1 

The mechanistic stuciy of PL binding provided information of 
varying molecular interactions with respect to decreasing distance 
from binding site (Figure 3A-E). The analysis of docking of AKTl 
has revealed that PL binds deep inside the cavity and is stabilized 
by the hydrophobic interactions (Figure 3F-G). However, the 
binding energy and dissociation constant of PL was lesser than the 
bound inhibitor (Table 2, Table 3). This may occur due to the 
smaller size of PL than the bound inhibitor. The Dock score was 
also highly negative and number of hydrophobic interactions that 
kept PL bound in the cavity was also reasonably high (27 
interactions from 6 different residues) as mentioned in Table 2. 
The key residue in AKTl for PL binding was Trp-80, which was 
common through all phases of PL (un)binding simulation, 
demonstrating its importance in initial binding of PL and finally 



bringing it into the active site (Figure 3). The Trp-80 was also 
involved in the majority of hydrophobic interactions and showed 
highest decrease in its solvent accessibility after PL binding 
(approx 60 A ) as shown in Table 3. This finding was consistent 
with a study [83] in which it was shown that the inhibition by Akti 
(an inhibitor of AKT) is critically dependent upon a solvent- 
exposed tryptophan residue (Trp-80) present in all three AKT 
isoforms and whose mutation to an alanine yields an Akti-resistant 
kinase. It was observed that the Trp-80 acts as the key interacting 
residue in docking and (un)binding simulation study of PL 
suggesting that PL also inhibits AKTl by the mechanism similar 
to that of Akti inhibitor and, therefore, validated our docking 
predictions for AKTl as well. The Akti acts as non-ATP- 
competitive allosteric inhibitor and blocks AKT phosphorylation 
and does not lead to active phosphorylated state of AKT [83] and 
thus the PL also inhibits AKTl by the similar mechanism i.e., acts 




Figure 4. Plumbagin (PL) binding to the cavity of Bcl-2. Panels A-D: show the (un)blnding simulation phases of PL, 'A' is farthest from the 
binding site, 'C is the closest to the binding site, and 'D' is the binding site phase. The residues involved in hydrophobic interactions are shown as red 
arcs. Those residues which are common to the last phase (D) are encircled. Panel E: Another representation for phase D. The whole protein is 
displayed in cartoon representation and the ligand molecules are in sticks; PL colored as green and the bound known inhibitor in blue. The 
interacting residues are labeled and shown as surface in different colors. 
doi:1 0.1 371/journal.pone.0087309.g004 
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Figure 5. Plumbagin (PL) binding to the cavity of NF-kB. Panels A-B: show the (un)blnding simulation phases of PL, 'A' is farthest from the 
binding site and representing many stages always coming with 2 interacting residues and 'B' is the binding site phase. The hydrogen bonds are 
shown as green-dashed lines with indicated bond length and the residues involved in hydrophobic interactions are shown as red arcs. Those residues 
which are common to the last phase (B) are encircled. Panel C: Another representation for phase B. The whole protein is displayed in cartoon 
representation and PL in sticks colored as green. The interacting residues are labeled and shown as surface in different colors. 
doi:1 0.1 371 /journal.pone.0087309.g005 



as non-ATP-competitive allosteric inhibitor. The inhibitory 
mechanism of PL (similar to Akti) is in contrast to ATP- 
competitive inhibitors, which protect the phosphorylated site of 
AKT from dephosphoryation by phosphatases [83] . 

Plumbagin docking and (un)binding simulation study of 
Bcl-2 

The (un)binding simulation phases of PL with increasing 
number of molecular interactions while approaching towards the 
binding site are shown in Figure 4A— D. The docking phase (Figure 
4D-E) showed molecular interactions of Bcl-2 with PL which binds 
on the exterior region of Bcl-2 cavity using hydrophobic 
interactions. Owing to the large size of previously known Bcl-2 
inhibitor (acyl-sulfonamide-based ligand) [84] in the complex 
structure (PDB ID: 2021), the binding constant and binding 
energy of PL was much smaller than that of the inhibitor (Table 2). 
In previous studies of molecular docking of inhibitor compounds, 
the residue Tyr- 199 of Bcl-2 was involved in hydrogen bonding 
[85,86]. However, in the present study, Tyr- 199 was not involved 
in hydrogen bonding, rather this residue was involved in 
maximum number of hydrophobic interactions and lost maximum 
surface accessibility (47 A^) while switching to bound state (Figure 
4D, Table 3). The Tyr- 199 was also found to be common in all 
phases of (un)binding simulation of PL which again shows its 
importance for PL binding (Figure 4). All the observations for Tyr- 
199, such as the maximum number of hydrophobic interactions 
involved, the maximum loss in surface accessibility, and the 
common residues in (un)binding simulation were corroborating 
with one another. Furthermore, the key residue involved, Tyr- 199, 
was also same as observed in the previous docking studies with 
other inhibitors [85-86]. 

Plumbagin docking and (un)binding simulation study of 
NF-kB 

The molecular docking and (un)binding simulation study of NF- 
kB with the interacting residues and molecular interactions are 
shown in Figure 5A-B. The first phase (Figure 5A) represents 



multiple stages of binding where the interacting residues were 
same but PL is approaching towards binding site and at variable 
distance. In the molecular docking phase (Figure 5B-C), PL 
packed against the residues Lys-123, Lys-122, Val-121, Tyr-36, 
and Asn-155. As shown in the Figure 5B, PL formed two hydrogen 
bonds, one between O-atom of 1 -carbonyl group of PL and N- 
atom of amide group of Asn-155 (2.98 A) and the other between 
O-atom of 4-carbonyl group of PL and ot-amino of Lys-122 (2.98 
A). There were 1 1 hydrophobic interactions from 4 different 
residues (Table 2, Table 3). It has been reported [87] that the Lys- 
122 and Lys-123 of NF-kB p65 subunit are the only acetyl- 
acceptor lysines which contact with DNA in minor groove as 
revealed from the crystal structure of the p50/p65 heterodimer 
bound to DNA. The acetylation of these residues reduces kB-DNA 
binding (may be neutralized by the positive charge on Lys-122 and 
Lys-123 s-amino groups) and represses its transcriptional activity, 
and thus leads to attenuation of p65-mediated transcription [87]. 
In our study, of the two lysine residues, Lys-122 formed H-bond 
and Lys-123 was involved in non-bonding contact with PL. Lys- 
123 was also found common among all phases of (un)binding 
simulation study (Figure 5). Therefore, the lysine residues (Lys- 
122, Lys-123) were not available for DNA binding, which resulted 
in the reduction in the DNA binding activity. The Dock score, the 
binding energy, and the dissociation constant were also compa- 
rable with the docking results of other chosen targets (Table 2). 

Plumbagin docking and (un)binding simulation study of 
Stats 

The (un)binding simulation phases with the increasing molec- 
ular interactions are shown in Figure 6A-C. The docking phase 
(Figure 6C-D) showed that PL bound to the DNA binding domain 
and packed against the residues Gly-421, Gly-419, Cys-418, Arg- 
417, Glu-415, and Arg-382. The Arg-417 was involved in 
hydrogen bonding interaction through the N^ atom of its 
guanidium group (Figure 6C) and also showed maximum loss of 
solvent accessibility in PL bound state (Table 3). Other residues 
involving hydrophobic contacts also showed significant loss in 
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Figure 6. Plumbagin (PL) binding to the cavity of StatS. Panels A-C: show the (un)blnding simulation phases of PL, 'A' is farthest from the 
binding site, 'B' is the closest to the binding site, and 'C is the binding site phase. The hydrogen bonds are shown as green-dashed lines with 
indicated bond length and the residues involved in hydrophobic interactions are shown as red arcs. Those residues which are common to the last 
phase (C) are encircled. Panel D: Another representation for phase C. The whole protein is displayed in cartoon representation and PL in sticks colored 
as green. The interacting residues are labeled and shown as surface in different colors. 
doi:1 0.1 371 /journal.pone.0087309.g006 

solvent accessibility in switching from unbound to bound state activity of StatS. The NHl of guanidium group of Arg-417 is 
(Table 3). A mutational analysis study [88] created the dual involved in hydrogen bonding with phospho-diester bonds of bases 
mutation Arg-4 14/4 17 that led to complete loss of DNA binding T-1007 and T-2007 (from complementary strand) [89-90]. 
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Similarly, another residue Arg-382, is also reported to be invok ed 
in hydrogen bonding with phospho-diester bonds of bases T-1006 
and T-2006 (from complementary strand) [89-90]. In another 
study, it has been reported that the residue Arg-382 when mutated 
to Trp and Gin abrogated the DNA-binding ability of Stat3 homo- 
dimer [91]. In our study, Arg-417 with 2 other residues (Gys-418, 
Gly-419) was found to be common throughout all phases of 
(un)binding simulation study (Figure 6). The residues Arg-417 (H- 
bonding, maximum loss of solvent accessibility, common in 
(un)binding simulation) and Arg-382 (hydrophobic contacts) play 
important role in PL binding, which is in agreement with the 
aforementioned studies. The Dock score is also highly negative 
comparable to the docking results of other chosen targets. The 
binding energy and the dissociation constant for PL were also 
similar to the results of other chosen targets (Table 2). 

Conclusions 

The present study for the first time used docking and 
(un)binding simulation analysis to identify intc-racting residues of 
five important cancer-signaling proteins and their molecular 
interactions with PL. Plumbagin inhibited the molecules by 
involving their functionally important residues leading to the loss 
of function. Among various important residues identified during 
the docking analysis, the key residues for the chosen proteins that 
play major role in PL mediated inhibition are: the PKKy residue, 
Val-882, involved in the hydrogen-bonding interaction; the AKTl 
residue, Trp-80, showing highest decrease in the solvent accessi- 
bility after PL binding; the Bcl-2 residue, Tyr-199, showing 
maximum hydrophobic interactions; the NF-kB residue, Lys-122, 
involved in H-bond formation; and the Stat3 residue, Arg-417, 
forming H-bond. The exact mode of inhibition was indicated 
when multiple modes of inhibition existed (i.e. AKTl). In all the 
five signaling target molecules, the Dock score was highly negative, 
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which showed good quality docking. We also presented the exact 
binding of PL and identified and characterized various interacting 
residues within the active sites of respective chosen proteins. In 

addition, we ranked the importance of residues for PL binding 
using ASA analysis. Also, the importance of residues was shown by 
ASA analysis (higher loss in solvent accessibility) and interaction 
analysis (high number of hydrophobic interactions or presence of 
H-bond) and by their consistent appearance in (un)binding 
simulation phases. A peculiar observation in this study was that 
PL formed hydrogen bond interactions with the target proteins 
only through the oxygen-atom of 1-carbonyl group which acted as 
a strong hydrogen donor. This finding may be useful for increasing 
drug potency in lead optimization where PL will serve as starting 
lead compound. The current study of PL docking analysis provides 
the structural insights into the binding mechanism of PL to the five 
chosen cancer signaling molecules. The PL-protein architecture is 
expected to guide scientists as a suitable model for understanding 
the multi-targeting drug structure and in providing structural 
details for the inhibitory mechanisms. 
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